Entangled Qubits in a non- Gaussian Quantum State 



T. Kiesel, W. Vogel 

Arbeitsgruppe Quantenoptik, Institut fiir Physik, Universitdt Rostock, D-18051 Rostock, Germany 

B. Hage 

ARC Centre of Excellence for Quantum-Atom Optics, Department of Quantum Science, 
The Australian National University, Canberra, Australian Capital Territory 0200, Australia 

R. Schnabel 

Institut fiir Gravitationsphysik, 
Leibniz Universitdt Hannover and Max-Planck-Institut fiir Gravitationsphysik (Albert-Einstein-Institute), 
Callinstrasse 38, 30167 Hannover, Germany 

We experimentally generate and tomographically characterize a mixed, genuinely non-Gaussian bipartite 
continuous-variable entangled state. By testing entanglement in 2 x 2-dimensional two-qubit subspaces, en- 
tangled qubits are localized within the density matrix, which, firstly, proves the distillability of the state and, 
secondly, is useful to estimate the efficiency and test the applicability of distillation protocols. In our exam- 
ple, the entangled qubits are arranged in the density matrix in an asymmetric way, i.e. entanglement is found 
between diverse qubits composed of different photon number states, although the entangled state is symmetric 
under exchanging the modes. 
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I. INTRODUCTION. 



Since the early days of quantum mechanics, entanglement 
has been a topic of great interest flT] |2l. Today, it is con- 
sidered as the key resource for many applications, such as 
quantum computation, quantum communication and quantum 
cryptography, for recent reviews cf. EllH. Last but not least, 
nowadays the relevance of entanglement is even discussed 
in the context of life sciences, cf. e.g. JSj. Traditionally 
there is a distinct treatment of entanglement in two different 
regimes, the discrete variable regime being based for instance 
on qubits, and the continuous variable (CV) regime having ei- 
ther Gaussian or non-Gaussian statistics, respectively. How- 
ever, both regimes are connected to each other, since qubits 
can be constructed out of CV states |6|, and vice versa entan- 
gled CV states may be composed of entangled qubits. 

A quantum state is refeiTed to as being entangled if it can- 
not be written as a statistical mixture of factorizable (uncor- 
related) quantum states Q- Although this definition is intu- 
itively clear, it is hard to check in practice for general quantum 
states. For the case of CV states the Peres condition verifies an 
important class of entanglement by partial transposition (PT) 
of the density matrix |8|. The special case of Gaussian entan- 
gled states can be fully characterized by moments of second 
order |9, 10|. Gaussian states are easy to prepare and it is 
easy to identify their entanglement. However, their applica- 
tion in quantum technology is limited. Non-Gaussianity of 
CV entangled states can be considered as a resource, which 
is of great interest for various applications in quantum infor- 
mation technologies. A general reformulation of the Peres 
criterion yields a hierarchy of inequalities in terms of observ- 
able moments of arbitrarily high orders 1 11 , 12 1. Recently, on 
this basis genuine non-Gaussian entanglement being invisible 
in second moments of the field quadratures, could be experi- 



mentally demonstrated fT3l . 

PT entanglement tests are necessary and sufficient for bi- 
partite entanglement in very special cases only: for Hilbert 
spaces of dimension 2x2 and 2x3 and for Gaussian C V quan- 
tum states, cf. e.g. EJH]. Entanglement of general states can 
be verified by involved generalizations of the PT tests lfT4l . 
or, based on entanglement witnesses[15|, by optimized entan- 
glement conditions 1 16 1. However, it has been proven that for 
any entangled state there exist subspaces of finite dimension 
in which the entanglement already exists |17|. A systematic 
search for these subspaces provides the ability to uncover the 
structure of the entanglement. In particular, the identification 
of these entangled subsystems may be helpful for designing 
protocols for the entanglement distillation, which extracts, in 
the ideal case, a maximally entangled state from the given 
mixed state, providing advantages in many quantum informa- 
tion tasks lH]. 

In the present article we study the entanglement of a 
CV quantum state, experimentally created by mixing a fully 
phase-randomized squeezed vacuum state by a balanced beam 
splitter with vacuum noise. We analyze the entanglement 
structure of this decohered nonclassical state by PT tests on 
two-qubit subsystems. In this way, we extend the complete- 
ness of the Peres criterion for bipartite qubit systems to the de- 
tection of entanglement in arbitrary continuous-variable sys- 
tems, and establish a link between entanglement in discrete 
variable and continuous variable systems. Several entangled 
qubits within the density matrix are identified with an over- 
whelming statistical significance. Surprisingly, the entangle- 
ment is only found between qubits of diverse nature. A closer 
inspection shows that the entanglement of these qubits can be 
distilled by the protocol of Bennett et al (BBPSSW) 1 18 1. 

The paper is structured as follows: In Section|llj we present 
the experimental setup and describe the state whose entangle- 



ment is analyzed. In Section III we examine the detection 
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of entanglement with Gaussian criteria, showing that they are 
not appUcable for our states. The search for entangled qubit 
subsystems is given in Section |IV[ and its use for choosing 
a suitable distillation protocol is discussed. A summary and 
some conclusions are given in SecfV] 



II. GENERATION OF A NON-GAUSSIAN ENTANGLED 
STATE. 

We start with a Gaussian squeezed- vacuum state being fully 
characterized by its variances Vx and Vp of the quadratures x 
and p, respectively, together with its orientation angle ip in 
phase space. To obtain a phase-independent non-Gaussian 
state, the orientation is uniformly distributed over a 2tt in- 
terval. This state is not squeezed anymore, but still exhibits 
significantly nonclassical properties 1 19 1, which are necessary 
for the generation of entanglement. In the second step, we 
send the phase-randomized state to a balanced beam splitter 
where it is mixed with a vacuum mode. For a sketch of the 
experimental setup, see Fig.[T] The two output modes, which 
form the entangled state examined in this article, are mea- 
sured by joint balanced homodyne detectors. This scenario 
is equivalent to propagate an entangled two-mode squeezed 
vacuum through a medium preserving the phase difference 
but destroying the absolute phase, similar to the experiment 
inRef. 120!. 

The squeezed mode was generated by a degenerate op- 
tical parametric amplifier (OPA). The OPA consisted of a 
type-I non-critically phase matched second order nonlinear 
crystal (7% Mg:LiNb03) inside a standing wave optical res- 
onator with a line width of 25 MHz. The OPA process was 
continuously pumped with second harmonic light at 532 nm. 
The resonator tuning was controlled via the common Pound- 
Drever-Hall method using phase modulated fundamental light 
at 1064 nm. For technical reason the pump phase was con- 
trolled such that the fundamental control field was deampli- 
fied, i.e. the OPA generated amplitude quadrature squeezing 
with respect to this field. With this setup we directly mea- 
sured a squeezing variance of -4.8 dB and an anti-squeezing 
variance of H-9.0dB, both with respect to the unity vacuum 
noise variance. 

The squeezed field propagated in free space from the OPA 
passing high-reflection mirrors. In order to apply the phase 
diffusion two of these mirrors were moved by piezo-electric 
transducers (PZT) driven by a quasi-random voltage. The 
driving signal was generated by a high quality PC sound card 
connected to appropriate filters and amplifiers. The sound 
card played back a previously generated sound file which was 
carefully designed to yield the desired phase distribution. 

Balanced homodyne detection (BHD) was used to measure 
the quadrature amplitudes of the two-mode state under con- 
sideration. The visibility with the spatially filtered local os- 
cillators were in the range of 98.5% to 99%. Regarding the 
detection phases the only meaningful figure is the difference 
between the detection phases of BHDa and BHDb because the 
state before the beamsplitter (BS) is phase randomized. The 
control of this phase difference was achieved by injecting an 




Figure 1. Simplified sketch of the experimental setup. LO: local os- 
cillator, OPA: optical parametric amplifier (squeezed light source), 
BHD: balanced homodyne detector, PZT: piezo-electrically actuated 
mirror applying the phase noise. Aux: Frequency shifted weak field 
providing a read-out for the difference of the detection phases of 
BHDa and BHDb. 



auxiliary field into the open port of the beamsplitter, which 
was frequency shifted with respect to the fundamental (LO) 
frequency. The demodulation of the BHD signal at the beat 
frequency provided an error signal for locking the detected 
quadrature angle to the phase of the auxiliary field with an off- 
set given by the electronic demodulation phase. The signals 
of the two individual BHD-photodetectors were electronically 
mixed down at 7 MHz and low pass filtered with a bandwidth 
of 400 kHz to address a modulation mode of the light showing 
good squeezing and a high dark noise clearance of the order 
of 20 dB. The resulting signals were fed into a PC based data 
acquisition system and sampled with one million samples per 
second and 14 bit resolution. For a more detailed description 
of the main parts of the setup we refer to 1211 l22ll . 

We obtained a set of 10^ quadrature pairs per detection 
phase configuration, for 10 equally spaced phases per output 
mode A and B. From these measurements, we estimated the 
density matrix elements of the state and their full covariance 
matrix via appropriate pattern functions, see Appendix [A| for 
details. Due to the central limit theorem, this yields the full 
statistical information about the experimental result. 

The reconstructed density matrix of the measured state is 
shown in Fig.[2]for the first 5^ density matrix elements in the 
Fock basis. Our reconstruction is restricted to the 7 lowest 
Fock states for each mode, thus the full bipartite density ma- 
trix has a rank of 7 x 7 = 49. The trace of the reconstructed 
state equals to 0.9836 ± 0.0002, such that the main informa- 
tion about the state is covered. The standard deviation of the 
single matrix entries is bounded by 0.0003. It is noteworthy 
that a higher-dimensional state reconstruction makes no sense 
for the entanglement test. Already a reconstructed matrix of 
rank 8x8 would consist of several entries that are no longer 
significant relative to the statistical error, as it is obtained for 
our sample. The result of the state reconstruction shows that 
all significant matrix elements are located on certain lines. 
Only matrix elements with k + m = I + n contribute, having 
their origin in the action of the beam splitter on the Fock state 
\k + m) in the input state. The phase diffusion eliminates all 
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Figure 2. Reconstructed two-mode density matrix of the examined 
entangled state in Fock basis. A block mn provides tlie absolute 
values of the density matrix elements pki,mn where k, I refers to the 
photon numbers at BHDa, and m, n to those at BHDb. 



coherences between Fock states in the input state. 



III. GAUSSIAN ENTANGLEMENT CRITERIA. 

Entanglement of Gaussian states can be completely char- 
acterized by Simon's inequality [10|. Here we extend the no- 
tion of Gaussian entanglement to all quantum states for which 
the entanglement can be identified by the condition of Simon. 
This can be also the case for non-Gaussian states. Hence en- 
tangled states which cannot be identified by this condition will 
be denoted as genuinely non-Gaussian ones. 

In the following we prove that Simon's condition for entan- 
glement is equivalent to squeezing of the input state. Let us 
start from a quadrature covariance matrix of a bipartite state 
in block form 



V = 



A C 
B 



(1) 



and define J 



The state itself does not have to 



with B being the covariance matrices of the subsystems, 
1 

-1 

be a Gaussian one, but we only examine its covariance matrix. 
A state is entangled if the following inequality is violated: 

det A det B+{\ - \ det C\f - TrAJCJBJC^J 

> lidet A + dct B) . (2) 

Let us assume that the covariance matrix of the initial state for 
the quadratures iin , Pin is given by 



(3) 



A beam splitter recombines the quadratures of this input field 
with the quadratures of vacuum to the quadratures of the two- 
mode output field = tXin + ^2:vac, 2^4 = — ''^in + i^vac, 



where the field transmissivity t and reflectivity r satisfy |ip + 
|rp = 1. Simon's criterion for the resulting covariance matrix 
leads to the following inequality for entanglement: 



t'{l-f)[iVp-l)iV,-'^)-C^ <0 



(4) 



Since < t'^ < 1, the term in square brackets has to be nega- 
tive. In the following we show that this condition is equivalent 
to squeezing of the input state. The eigenvalues of the co- 
variance matrix of the latter are the roots of the characteristic 
polynomial. 



pix) = {v, - x){Vp - X) - a 



2 

xp' 



(5) 



The two roots Ai_2 are the minimum and maximum quadrature 
variance of the state. A state is squeezed if for one of these 
roots we find Ai < i, while for the other we have X2 > ^■ 
This is the case if and only if p{X) is negative between both 
roots, i.e. p{-k) < 0. This leads to 



p(i) = (V^,-i)(Fp-i)-C2 <0, 



(6) 



which is equivalent to Simon's condition for entanglement. 
Therefore, a beam splitter creates Gaussian entanglement if 
and only if the input state is squeezed. In this sense, one may 
state that Gaussian entanglement of the split state has its origin 
in Gaussian nonclassicality of the input state. We emphasize 
that this fact is not restricted to Gaussian states, but holds for 
arbitrary quantum states. 

One can proceed with the analysis of higher-order mo- 
ments. In fT3l, a criterion containing fourth-order moments 
was sufficient to demonstrate genuine non-Gaussian entangle- 
ment. For the studied states the Simon test, based on second 
order-moments, failed. In our case, one can show that one has 
to go to sixth-order moments for this purpose, which is a cum- 
bersome task. Therefore, we propose an alternative way to 
verify non-Gaussian entanglement, which does not only ver- 
ify the entanglement. In addition it provides useful insight 
into the entanglement structure of the state. 



IV. ENTANGLED QUBITS. 

Recently it has been been proven that any entangled quan- 
tum state must also be entangled in a finite dimensional 
Hilbert space |17|. Therefore, it is sufficient to detect entan- 
glement in a submatrix of the full density matrix. The Peres 
criterion states that for any separable state the partial transpose 
of the density matrix, defined by its matrix elements 



{k, m\ p \l, n) — (fc, n\ p \l, m) 



(7) 



is positive semidefinite fSl. Hence, we may start to project 
the state onto an arbitrary two-qubit subsystem and calculate 
the eigenvalues of the partially transposed density matrix. If 
at least one of the latter is negative, the state has a negative 
partial transposition and entanglement has been verified by 
the Peres-criterion. This test can be implemented efficiently, 
since a bipartite qubit subsystem can be described by a 4 x 4 
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Figure 3. Probabilities for the occurrence of entangled two-qubit sub- 
spaces. The statistical significance of the successful entanglement 
tests is color-coded. It represents the smallest (negative) eigenvalue 
of the partially transposed state in units of the corresponding standard 
deviation a. 



matrix, which has only four eigenvalues. The statistical un- 
certainty (7(Amin) of the eigenvalue Amin can be obtained by 
a Monte-Carlo simulation: We draw random matrices, whose 
entries are chosen by a Gaussian distribution around the re- 
constructed state, with the covariance matrix as determined 
from the experimental data. With these simulated density ma- 
trices, which are consistent with the reconstructed state within 
its statistical uncertainties, we calculate the set of least eigen- 
values. From the statistics of the results we can estimate the 
standard deviation of our experimental entanglement test. We 
note that whenever such a test fails for two-qubit subsystems 
it may be successful in subsystems of higher dimensions. 

This test has been performed for all possible qubit subsys- 
tems. In Fig. [3] all two-qubit subsystems with negative partial 
transpose are illustrated by the probability of their occurrence 
in the state and by the statistical significance, ^ , of the 

smallest eigenvalue. For each mode, the qubit subsystems are 
ordered from |0) , |1) to |0) , |6), continuing with |1) , |2) to 
|1) , |6) and so on. Each of the two axes running over sub- 
system labels refer to one of the two modes. First, it is obvi- 
ous that there is no symmetric subsystem which displays the 
entanglement. That is, if both modes are projected onto the 
same subsystem, no entanglement is indicated. However, pro- 
jecting on different subsystems, such as the ones composed of 
the states |0) , |1) in one mode and |1) , |2) in the other, gives 
states whose entanglement can be verified with significance 
up to 128 standard deviations. In total we find 10 asymmet- 
ric pairs of qubits whose entanglement has a statistical sig- 
nificance of more than two standard deviations. The asym- 
metry is remarkable, since the state itself is symmetric with 
respect to both modes. The knowledge of the structure of the 
entanglement of a given mixed quantum state, in particular the 
identification of the localization of entangled qubits within a 
complex CV state, is essential for applications. 

Since we find entangled two-qubit subsystems in the state 
under study, our results clearly show that entanglement dis- 
tillation is possible. Now the question for a suitable proto- 



Figure 4. Probability of finding a singlet state \tj) ) in a specific 
qubit subsystem. 



col arises. One scheme, which works for any two-qubit state, 
has been given in [23 1. Our method as presented here identi- 
fies appropriate two-qubit Hilbert subspaces. Any subspace in 
which entanglement is found, may be chosen for distillation. 
The higher the probability of occurrence of the subsystem in 
the full state, the more efficient is the distillation. 

For some entangled qubits, one can directly apply the 
BBPSSW protocol [ISJ- Let us have a look at the fidelity of a 
singlet state with respect to the experimentally generated one 
within the chosen two-qubit subspace. A singlet in this two- 
qubit space is defined as 1^") = ^ (|iA,Js> - |jA,»s>), 
where , and jis) , |j\b) are Fock states forming a ba- 
sis for the qubit systems in the modes A and B, respectively. 
If the fidelity between the state and the singlet, 

F={^-\p\^-), (8) 

is greater than the BBPSSW protocol can be directly ap- 
plied to distill the state. In Fig. [4j the fidelities for the subsys- 
tems of interest are shown. Obviously, several subspaces have 
a fidelity which is significantly larger than i being suitable for 
the application of the BBPSSW protocol. 



V. SUMMARY AND CONCLUSIONS 

We have demonstrated an efficient method to examine the 
entanglement structure of CV quantum states using the exam- 
ple of a mixed, genuine non-Gaussian entangled state whose 
entanglement is invisible for Gaussian entanglement criteria. 
We have analyzed the two-qubit subsystems and found neg- 
ative partial transpositions with large statistical significances 
of up to 128 standard deviations. Our method proves that the 
CV state under consideration is entangled and distillable. 

A remarkable result of our analysis is that the decohered 
state under consideration only shows entanglement in two- 
qubit subsystems composed of diverse qubits, for instance 
|0) , |1) and |2) , |1). Such an insight into the entanglement 
structure might be used to analyze the decoherence process 
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that turned an initially pure entangled state into a mixed en- 
tangled one. Moreover, the knowledge of the entanglement 
structure is important for the development of proper filtering 
methods for the aim of entanglement distillation. 
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~{r,i) 
PkLmn 



Plnn]klij^U,jB) 



(A3) 



Appendix A: Reconstruction of the density matrix. 



where the brackets are symbols for 



We obtained a set of N = 10^ quadrature pairs 
{xj{4'j^),yj{(j)ji^)}jLi data points per phase configuration, 

for = 10 equally spaced phases (pj^ ^ = '^^^^ per out- 
put mode A, B. From these measurements, we estimated the 
density matrix of the state via appropriate pattern functions, 

xe^('"-")*^«/™„(y,(0,J)},(Al) 
xe^('"-")*^-/™„(%(0,J)},(A2) 



N 



{Fi3,jA,jB)) 



1 

NjN 



3A,JB) 



(A4) 



Of course, all diagonal elements of the density matrix are real: 

(i) 

Pkk mm = 0. Furthermore, the covariance matrix of all entries 
can be estimated in the standard way as 



'^^^ \ Pkl,mnPk'l',n 



1 



NiN 



F 



rnn.kl 



( j, JA , jB)Fl^:n',k'l' (i' JA , 3B ) 



Pki Pk 'V ,m'n' 



(A5) 



Since the quadrature pairs are identically and independently 
distributed for any fixed pair of phases 0j ^ , (pj^ , the sums 
over the pattern functions are Gaussian distributed due to the 



central limit theorem. Therefore, Eqs ( A3 1 and ( A5 1 provide 
the full statistical information about the sampled density ma- 
trix. 
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